The nonlocal correlation: A key to the solution of the CO adsorption puzzle 
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We study the chemisorption of CO molecule into sites of different coordination on (111) surfaces 
of late 4d and 5d transition metals. In an attempt to solve the well-known CO adsorption puzzle 
we have applied the relatively new vdW-DF theory of nonlocal correlation. The application of the 
vdW-DF functional in all considered cases improves or completely solves the discrepancies of the 
adsorption site preference and improves the value of the adsorption energy. By introducing a cutoff 
distance for nonlocal interaction we pinpoint the length scale at which the correlation plays a major 
role in the systems considered. 

PACS numbers: 68.43.Bc, 68.47.Dc, 71. 15. Mb, 71.45.Gm, 31.15.cs,34.35.+a 



For decades the density functional theory (DFT) [l| 
has been the main theoretical tool used to analyze, un- 
derstand and predict material properties and chemical 
processes. Although very successful in explaining the 
finest details of many systems, even very complex ones, 
it fails in some cases. In particular it fails in weakly 
bound systems where van der Waals forces play a major 
role. Examples are numerous. Practically all biologi- 
cal systems and any system containing noble gas atoms 
cannot be calculated correctly within the current approx- 
imations to the DFT, the local (LDA) and the semilocal 
generalized gradient approximation (GGA). A practical 
solution to this kind of problems was relatively recently 
proposed by Dion et al. Q ■ They have developed a new 
exchange-correlation functional named vdW-DF, which 
is parameter free and therefore is considered ah initio. 
The main idea of the vdW-DF functional is the intro- 
duction of nonlocal correlation into DFT calculations. 
The failure of local and semilocal functionals in weekly 
bounded systems is expected, but what is puzzling is that 
current approximations to the DFT fail to predict the 
correct adsorption site of CO molecules on several (111) 
metal surfaces, of which the Pt is the most famous ex- 
ample. Adsorption of CO on Pt surface is of tremen- 
dous practical importance, in particular in exhaust cat- 
alyst in the car industry. Some eight years ago Feibel- 
man et al. |j| published a comprehensive study of CO on 
Pt(lll), performing calculations independently with sev- 
eral state-of-the-art DFT computer codes, for a variety 
of exchange-correlation potentials. They found that DFT 
systematically predicts that CO prefers the hollow site on 
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Pt(lll) surface, while experiments show the adsorption 
into the top, sometimes bridge but never hollow site. Due 
to expected success of DFT in such chemisorbed systems 
the problem was named the CO adsorption puzzle. The 
same puzzle appears on some other (111) surfaces such 
as Cu, Rh, Ag, and Au. Even more confusing is the cor- 
rect prediction of the hollow site on some other surfaces 
such as Pd(lll) in agreement with experiment. Addi- 
tionally to the correct site preference getting the correct 
barrier height between the sites of different coordination 
is of tremendous importance for the analysis of diffusion 
phenomena on surfaces as discussed by Hu et al. Q . By 
considering the puzzle problem they have also noticed 
that the lack of nonlocal correlation is one of the possi- 
ble sources of the problem. 

In order to solve the puzzle we apply the vdW-DF 
functional. Even if up to now vdW-DF was mostly used 
to treat van der Waals type of systems, achieving quite 
good success [I| , the functional itself was constructed in a 
seamless fashion, i.e. describing the nonlocal correlation 
correctly at all length scales. In fact, the functional has 
recently been successfully applied to a covalently bonded 
system @. We apply the vdW-DF functional as a post- 
processing procedure on a charge density obtained from 
a standard GGA calculation. Such an application of the 
functional has been justified by the work of Thonhauser 
et al. 7] , who implemented the functional in a fully self- 
consistent manner into a DFT code and showed that the 
change of Kohn-Sham eigenstates was negligible. 

Following the original procedure of Dion et al. we re- 
place the conventional GGA correlation energy by a sum 
of two terms, the LDA correlation and a newly defined 
nonlocal term E£ ■ 

-BvdW-DF = EGGA — -BGGA.c + ^LDA.c + (1) 

where Eqqa is the selfconsistent GGA total energy. The 
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last term in Eq. ([T]) is completely determined by the 
electronic charge density and requires the evaluation of 
a double space integral: 
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where n(r) is the charge density at point r and <fi(r, r') 
is a kernel function described in Ref. Q. Since to our 
knowledge there is no available code to calculate the 
vdW-DF functional we have developed our own fgj . The 
code has been tested on the few examples which have 
been published so far 0], achieving a complete agree- 
ment. Our ab initio DFT calculations were carried out 
mostly by employing the pseudopotcntial plane- wave for- 
mulation as implemented in PWscf 0. In the spirit of 
the original puzzle article Q we performed calculations 
for the experimentally most interesting coverage of 1/4 
CO monolayer (ML) using two other codes, the Dacapo 
program [10, llll, which is an ultrasoft pseudopotential 
plane- wave code similar to PWscf, and VASP [HEEEi, 
which is the most accurate of the three due to its PAW 
implementation (TEj . 

It has been argued by the authors of the vdW-DF func- 
tional and others 0, @| that the vdW-DF should be used 
together with the revPBE functional for the calculation 
of the exchange energy contribution, because it seems to 
approximate most correctly the exact exchange and by 
that avoids some spurious binding effects. We have used 
both PBE and revPBE [l]} functional in our DFT 
calculations. 

We have applied the vdW-DF functional to the ad- 
sorption of CO into top and fee hollow sites on a whole 
range of (111) surfaces which have been related to the 
CO puzzle, for two different coverages. The results are 
summarized in Table |TJ We immediately notice the im- 
portant result that in almost all cases the energy of ad- 
sorption into the top site is lowered by roughly 100 meV 
compared to that of the hollow site. In the case of Rh 
this completely solves the puzzle, while for Pd the hollow 
site remains more stable, which is in agreement with ex- 
periment. The shift in energy is somewhat smaller for Pt 
and the preference for the hollow site remains, confirming 
its reputation as the most difficult puzzle system. Fur- 
thermore, in all cases the adsorption energy is improved 
compared to the pure GGA value. 

In order to understand the difference between the PBE 
(or any GGA) and the vdW-DF correlation energy we 
have examined the spatial distribution of the correlation 
energy densities. For both we have calculated the distri- 
bution for the CO chemisorbed on metal, and subtracted 
the values for the clean surface and an isolated molecule 
(at exactly the same positions as in the coupled system), 
in order to find in which regions of space changes occur 
upon chemisorption. Since the PBE correlation depends 
upon the gradient of the charge density and the vdW- 
DF correlation upon the charge density itself through 



the integral ([2]), their spatial distribution is quite differ- 
ent. In Fig. Q] we show the difference of the two. In the 
red regions the PBE correlation overbinds compared to 
vdW-DF, while the opposite is true in the blue regions. 
Integration over space shows that the total effect is al- 
ways an overbinding by PBE, but it is more pronounced 
in the fee hollow case than in the top case, by about 
100 meV, as discussed in Table Q] The figure suggests 
that this is due to the fact that there are three regions 
of overbinding in the fee hollow case and only one in the 
top case (although the latter is somewhat stronger than 
any single hollow one). 




FIG. 1: Differences in correlation energy densities between 
revPBE (£ rcv pBE,c) and vdW-DF (E L oA,c + Ef) calculations, 
for top site, upper panel, and hollow site, lower panel. The 
three-dimensional representation (left) and the cut along the 
CO molecule and through (the nearest) metal atom of the 
Pt surface (right) are shown. Red (blue) colors show regions 
in which PBE yields larger (smaller) binding energies com- 
pared to the vdW-DF functional. Color scales are the same 
in both images, with green representing zero value. By sum- 
ming all the contributions we found out that PBE correlation 
overbinds compared to vdW-DF in both cases. To get insight 
into this figure please follow the links: 



|http://www.youtube.com/watch?v=IbPhmFb-Vsg 



|http://www.youtube.com/watch?v=rq xwZOMaWbg 



We conclude that the PBE correlation functional 
makes a quantitative error in evaluating correlation en- 
ergies of different adsorption sites in these systems and 
therefore seems to yield wrong results, e.g. due to the 
failure to correctly predict adsorption sites in seemingly 
chemically bonded systems — this was the puzzle. PBE 
otherwise correctly captures the physics of the problem 
as our relaxation tests for puzzle systems of Pt and Cu 
show. Since we do not have access to forces which would 
appear in a selfconsistcnt implementation of the vdW- 
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TABLE I: Comparison of the energies of CO adsorption on (111) surfaces of several metals, calculated with the standard 
GGA and the vdW-DF post-processing approach. The abbreviation in the first column denotes the program used: V-VASP, 
D-Dacapo, P-PWscf. The letter "r" refers to the selfconsistent use of the revPBE functional, while otherwise PBE (VASP, 
PWscf) or PW91 [ltj (Dacapo) were used. Columns 2-4: The energy differences between the top and the hollow site. The 
results marked by an asterisk are obtained using an exchange functional different from the one with which the electronic 
density was calculated. Column 5: The 'NL shift', i.e. the change between the values in column 2 and column 3 or 4 (whichever 
corresponds to the selfconsistent functional). Columns 6-8: The calculated adsorption energies for the top site, column 7: The 
experimental value, where available. In all cases CO molecule is at GGA equilibrium position. 



DF theory, we have relaxed the CO molecule by hand. 
We first changed the position of the molecule relative to 
the surface in the perpendicular direction while keeping 
the C-0 distance at the relaxed value of the GGA calcu- 
lation, and later we also changed that distance. It turns 
out that the position of CO determined in GGA calcula- 
tions is the optimal one also for the vdW-DF functional. 
Moreover, the curvature of the adsorption energy around 
the minimum does not change if wc apply vdW-DF the- 
ory although the value of the minimum shifts. In systems 
in which the PBE does not describe the physics correctly, 
as for example a Kr dimer, both parameters — the posi- 
tion of energy minimum and the curvature around it — are 
strongly changed upon the application of the vdW-DF 
theory [lo|. 

We have also made some calculations of CO adsorp- 
tion on (111) surfaces of noble metals Au and Ag, and 
found that the results are qualitatively different to the 
other considered surfaces. The experimental chemisorp- 
tion energies of CO on these surfaces are in the range 
of 300-400 meV. In comparison, pure GGA energies are 
somewhat too small, which is more reminiscent of true 
van der Waals systems where GGA clearly underbindes, 
than of the strongly bonded systems considered in this 
paper. Also, while for other systems the equilibrium po- 
sition remains the same after applying the vdW-DF func- 
tional, we found that for Au and Ag the position of the 
minimum moves outwards, and (in agreement with ex- 
periment) the top site becomes clearly preferred. These 
properties single out the noble metals as rather different 
from the other cases considered, and we leave the detailed 
study for the future. 



In order to further investigate how the nonlocality in 
the vdW-DF correlation compares with the usual semilo- 
cal GGA approach we have introduced a cutoff distance 
in our calculation of Eq. by a function 

e(d cut - \r - r'\) (3) 

under the integral. The cutoff d cut defines the largest 
allowed distance of the nonlocal correlation. 

We performed calculations for CO chemisorbed on 
Pt(lll) in top and hollow site at equilibrium distance, 
and subtracted the energy for the separated molecule 
and surface, all calculated with the same value of d cu t- 
The results are shown in Fig. [SJ Three regimes are im- 
mediately apparent. At small d C ut> less than ~ 1 A, 
the contribution to the nonlocal correlation comes from 
within regions where charge is accumulated or depleted 
upon chemisorption (roughly corresponding to regions of 
intense colors in Fig. [T]). This region is somewhat un- 
physical, because in the integral © only a small part of 
the kernel <f> contributes. In the region ~ 2 - 3 A. which 
is the characteristic length of the chemical bond (and of 
the distances between blue blobs in Fig. [J) the nonlocal 
correlation becomes more binding. These changes oc- 
cur faster for the hollow than for the top site, reflecting 
the fact that there are three bonding regions instead of 
one. Finally, in the true van der Waals region, say for 
d C ut > 4 A, the interaction is, of course, attractive, but 
this region contributes little to the strongly chemisorbed 
system of CO on Pt(lll), and almost nothing to the dif- 
ference between the top and hollow site. 

Fig- HI shows that the physical length scale of the non- 
local correlation energy, i.e. when the energies and the 
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FIG. 2: The energies of adsorption of a CO molecule at top 
and hollow sites on the Pt(lll) surface at a coverage of 1/12, 
and the difference of the two. The calculation was done using 
the density obtained in a revPBE GGA calculation, vary- 
ing the cut-off in Eq. <[3j . At d cut = (i.e. including GGA 
revPBE exchange but only the LDA term in correlation) the 
top site is clearly preferred. The full revPBE+NL (large d cu t) 
gives a small preference for the top site, while the plain GGA 
shows a strong preference for the hollow site, at variance with 
experiment. 



energy differences saturate, is of the order of the distance 
between the carbon and the substrate atom. However, at 
this length scale the exact geometry of the chemisorbed 
complex becomes relevant, which intuitively explains why 
the GGA, which depends only upon the local value of the 
density gradient and thus does not see the big picture, is 
not able to give the correct values, in particular the en- 
ergy differences between sites of different coordination. 

In order to understand more thoroughly how nonlocal 
correlation binding energy depends on d cu t we have vi- 
sualized the correlation energy density for two different 
values of parameter d cu t , for top and hollow site adsorp- 
tion, as shown in Fig. [3J 

From Fig. [3] we conclude that both intra- and inter- 
region contributions exist. Comparing Fig. [2] and Fig. [3J 
one can see that inter-region interaction starting to take 
place at d cu t ~ 2A is the most important nonlocal con- 
tribution to the energy difference between top and hol- 
low site, and that it is also responsible for the change 
of the adsorption energy. Looking carefully in the upper 
panel of Fig. [3Jit is noticeable that most of the nonlocal 
correlation binding energy in the top site adsorption is 
located just above the Pt metal atom (region 1), while in 
the hollow site case most of it is found just underneath 
the carbon atom (region 1), reflecting different chemical 
bonding mechanism in the two cases. 

The fact that the PBE functional as well as other GGA 
functionals ovcrbind the adsorption of molecules on metal 
surfaces was realized some time ago [2lj. In the whole 
letter we made little reference to the exchange energy, 




FIG. 3: (Color online) Binding energies due to NL term calcu- 
lated with infinite d cu t, upper panel, and with d cu t~l A lower 
panel. Left (right) column shows top (hollow) site adsorption. 
Color scales are the same in all four images, blue color repre- 
sents zero value. On the lower right image the sphere of diam- 
eter of 3 A is shown while the spheres representing Pt atoms 
have 1 A in diameter. From the upper images it is clear that 
the most of the binding energy is coming from two regions in 
the top and four regions in the hollow case which are numer- 
ated in the figure. What remains unclear from the upper panel 
is whether the correlation binding energy in marked regions 
is due to intra-region or inter-region interaction. (Compar- 
ing these images to the ones in Fig. [1] one can see that these 
regions correspond to the ones where PBE shows underbind- 
ing compared to vdW-DF, and that they are actually regions 
where the charge density was depleted due to binding). To an- 
swer this question we have performed calculations with cutoff 
parameter value d cu t = l A which is just large enough to take 
all intra-region interaction into account, while avoiding any 
inter-region interactions (for the length scale reference in the 
figure the spheres are drawn) . If there is no inter-region inter- 
actions the corresponding images from the upper and lower 
panel would not differ. However, a large difference between 
images exists, proving the presence of strong inter-region in- 
teractions. To get insight into this figure please follow the 

links: 

http://www.youtube.com/watch?v =kgMNEulddoo 
http://www.youtube.com/watch?v=6iuFD81hucY 



but followed the arguments of the authors of the vdW- 
DF theory, and considered the revPBE formulation as 
a possible improvement. However, we believe that the 
implementation of the exact exchange might be the best 
approach, and may even lead to the resolution of the 
puzzle. 

In conclusion, we have shown that the application of 
vdW-DF theory resolves a large part of the CO adsorp- 
tion puzzle. Even though in almost all cases we achieve 
a preference of the top site, we cannot claim a complete 
solution of the puzzle, because larger energy differences 
arc required to agree with experimental thermal stability. 
We have shown that nonlocal correlation effects at a typ- 
ical length scale of C-metal atom distance (~ 2 — 3 A) are 
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decisive for a good description of CO puzzle systems. We 
arc not sure whether the construction of a better semilo- 
cal functional is impossible, as was recently concluded on 
the grounds of exploration of the puzzle problem pjj , but 
in the construction of such a semilocal functional (which 
would be desirable due to a lower computational cost) 
one should take as a guide and the ultimate check of the 
quality of the functional the treatment of the same CO 
puzzle systems with vdW-DF, as it reveals the delicate 
balance between nonlocal correlation and adsorption site 
preference. To summarize, our calculations show that the 
vdW-DF theory is crucial in chemically bonded systems 
which were initially not expected to be much affected by 
it. 
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